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Abstract 

•rH , 

^ . We investigate the relationship of "physical" parton densities defined by fc-factorisation, to 

those in the minimal subtraction scheme, by comparing their small-x behaviour. We first 
summarize recent results on the above scheme change derived from the BFKL equation at 
NLx level, and we then propose a simple extension to the renormalisation-group improved 
(RGI) equation. In this way we are able the examine the difference between resummed 
gluon distributions in the Qo and MS schemes and also to show MS scheme resummed 
results for P gg and approximate ones for P qg . We find that, due to the stability of the 
RGI approach, small-x resummation effects are not much affected by the scheme-change 
in the gluon channel, while they are relatively more sensitive for the quark-gluon mixing. 



Predictions of perturbative Quantum Chromodynamics for DGLAP [1] evolution kernels in 
hard processes have been substantially improved in the past few years, both by higher order 
calculations for any Bjorken x [2] and by resummation methods in the small- a; region [3-14]. 
However, higher order splitting functions are factorisation-scheme dependent: while the NNLO 
results and standard parton densities [15, 16] are obtained in the minimal subtraction (MS) 
scheme, the resummed ones are in the so-called Qo-scheme, in which the gluon density is defined 
by fc-factorisation of a physical process. Therefore, in order to compare theoretical results, or 
to exploit the small- x results in the analysis of data, we need a precise understanding of the 
relationship of physical schemes based on fc-factorisation and of minimal subtraction ones, with 
sufficient accuracy. 

The starting point in this direction is the work of Catani, Hautmann and one of us (M.C.) [17, 
18], who calculated the leading- log a; (LLx) coefficient function R of the gluon density in the 
(dimensional) Qo-scheme 1 versus the minimal subtraction one, namely 
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where u is the moment index conjugated to x, the MS density 
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factorises a string of minimal-subtraction 1/e poles starting from an on-shell massless gluon, 
^(ds/uj) is the LLx BFKL [21] anomalous dimension, and the explicit form of R will be given 
shortly. 

The purpose of the present Letter is to show how to generalise the relation (0) to possibly 
resummed subleading-log levels and to quarks. We first summarize the essentials of such gen- 
eralisation at next-to-leading log(x) (NLx) level, following the detailed analysis of the BFKL 
equation in 4 + 2e dimensions of two of us [20]. We then propose a simple extension of the 
method to the renormalisation group improved (RGI) approach [4,6]. On this basis, we show 
the effect of a Q to MS scheme change on a toy resummed gluon distribution and we calculate 
a full MS small- x resummed P gg evolution kernel as well as a small-x resummed P qg evolution 
kernel in an approximation to the MS scheme. 



1 Scheme change to the MS gluon at NLcc level 

Let us first summarize the results of [20] for the gluon channel only (Nf = 0). The starting 
point is the BFKL equation [21] with NLx corrections [22, 23] continued to 4 + 2e dimensions, 
as described in more detail in [20] . In particular, we consider running coupling evolution at the 
level of the one-loop /3-function 

p(a s , e) = ea s - ba 2 s , b = (3) 

1Z7T 



so that 
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1 The label Qq referred originally [19] to the fact that the initial gluon, defined by fe-factorisation, was set 
off-mass-shell (k 2 — Qq) in order to cutoff the infrared singularities. It turns out [20], however, that the effective 
anomalous dimension at scale fc 2 3> Qq is independent of the cut-off procedure, whether of dimensional type or 
of off-mass-shell one. 
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where a M = (<7/x e ) 2 e _e ^ 1 )/(47r) 1+e is normalised according to the MS scheme. 

Note first that, the ultraviolet (UV) fixed point of Eq. Q at a s = e/b separates the evolu- 
tion into two distinct regimes, according to whether (i) a M < e/b or (ii) a^ > e/b. In the 
regime (i) a s (t) runs monotonically from a s = to a s = e/6 for — oo < t < +00 — and is thus 
infrared (IR) free and bounded — , while in the regime (ii) a s (t) starting from e/b in the UV 
limit, goes through the Landau pole at = log(l — e/ba^/e < 0, and reaches a s = from 
below in the IR limit. 

The main result of [20] is a factorisation formula for the BFKL gluon density in 4 + 2e 
dimensions. If the gluon is initially on-shell, and the ensuing IR singularities are regulated by 
e > in the (unphysical) running coupling regime (i) mentioned above, then the gluon density 
at scale k 2 = /i 2 e* factorises in the form 

gi Qo \t,uj) = K(a s (t),uj)exp [ dr y(a s (r), u; e) , (5) 

J — 00 

where 7 is the saddle point value of the anomalous dimension variable 7 conjugated to t and 
the factor 3\T e — which is perturbative in a s (t) and e — is due to fluctuations around the saddle 
point. The expression of 7 for e > is determined by the analogue of the BFKL eigenvalue 
function, namely at NLx level by the equation 



a 8 (t) 
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where, by definition, 



34°)(fc 2 r- 1 - e = xf ) (7)(fc 2 ) 7 - 1 , (7) 

3CW(fc 2 r 1 - a - = rf ) (7)(fc 2 ) 7 - 1 , (8) 

and the detailed form of the kernels is found in Refs. [22,23] on the basis of Refs. [24,25]. The 
result (0) is proven in [20] from the Fourier representation of the solutions of the NLx equation 
by using a saddle-point method in the limit of small e = O (ba s ), which singles out 7 as in 
Eq. ®. 

Let us now make the key observation that — due to the infinite IR evolution down to 
a s (— 00) = — the exponent in Eq. © develops 1/e singularities according to the identity 

/•««(*) da 

dr ^(a s (r),uj;e) = — — y(a,w,e) , (9) 

J a(e - ba) 

which produces single and higher order e-poles in the formal ba/e expansion of the denominator. 
However, such singularities are not yet in minimal subtraction form, because we can expand in 
the e variable the leading and NL parts of 7, as follows: 

7Kc;.)=7 (0) (^,e)+« s 7 (1) (^,e) (10) 
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While the e = part is already in minimal subtraction form, the terms O (e) and higher need 
to be expanded in the variable e — ba in order to cancel the series of e-poles generated by the 
denominator: 

e ba . 1 1 . 

+ 1, (11a) 

+ ba + e, (lib) 



e — ba e — ba 
e 2 b 2 a 2 
e — ba e — ba 



3 



and similarly for the higher order terms in e. Therefore, by replacing Eqs. (fTU|) and (fTT]) into 
Eq. (jSJ), we are able to factor out the minimal subtraction density in the form 
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where now the ^-independent MS anomalous dimension is 

7 (MS) (a s ,u;) = 7(a s ,u;;6a s ) 
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and contains some NNLx terms related to the e-dependent ones in square brackets in Eq. 



Correspondingly, the coefficient R £ in Eq. (|12jl has a finite e = limit, at fixed values of 



a,. 



and a s (t) = a M /(l + ba^t). Therefore, we are able to reach the physical UV-free regime (ii), 
and we obtain 
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which is the result for the R factor at NLx level we were looking for. 

A few remarks are in order. Firstly, the expansion coefficients 'j^ and 72°'* are simply 
obtained from the known form of the e-dependence of (7) = Xo{l) + £ Xi(l) +£ 2 X2{l) + O (e 3 ) 
in Eq. (JJJ), while 7^ is not explicitly known, because the e-dependence of x^e\l) in Eq. (JHJ) 
has yet to be extracted from the literature [24,25]. We quote the results 
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In particular 7!°'' , together with the LLx form of 

X ' 



7l V / -X' (0) (7l 



yields the result of Ref. [17, 18] 
Rv(a s (t),u)=R( lL (^l) 
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while 73 and 7^ provide the new NLx contribution to R of [20]. 
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Secondly, the anomalous dimension in the Qo-scheme takes contributions from No only, 
namely 2 

7 <« K a) = 7 (o» (I) + Q , 7 a) (|) _ _ (20) 

and is therefore independent of the kernel properties for £ > 0. On the other hand, by Eqs. (|13p 
and (UU), 7^ MS ^ is related to 9\ by the expression 



7 ^ K w) = T ^ + as7 u, ^ + to r , (2D 

whose origin is tied up to the identity (jlljl . Indeed, we have separated terms of order e or 
e 2 into minimal subtraction and coefficient contributions: therefore, their t-evolution should 
cancel out in the e = limit, which is the content of Eq. (j21j) . 

Using Eqs. (|20|) and (|21|) . the well known relations [17] for NLx anomalous dimensions can 
be extended to subleading levels, as generated by the ^-expansion. Thus the difference 



7 (MS)_ 7 (Qo) =6as 



.,(0) , (i) , , (o) . d log N 



7i + a s 7i + 6a s 72 



9 log a s 



(22) 



is computed up to NNLx level as outlined above, even if the dynamical e = NNLx contribu- 
tions to the 7's are not investigated here. 

We conclude that, while the anomalous dimension in the Qo-scheme (which is roughly a 
"maximal" subtraction one) only depends on the e = properties of the BFKL evolution, the 
MS coefficient and anomalous dimension both depend on higher orders in the e-expansion of 
the kernel eigenvalue, which generate subleading contributions. The result in Eq. (j22j) of [20] 
directly provides the scheme change for the gluon anomalous dimension at NNLx level. 



2 Resummed results for the MS gluon splitting function 

A problem exists concerning the magnitude of the scheme change summarized above in the 
small-x region. In fact, the explicit form of the coefficients X, 7{°' ) , 7^ , 72 exhibit leading 
Pomeron singularities of increasing weight, indicating that a small-x resummation is in principle 
required for the scheme-change too. As a consequence, any resummed evolution model should 
provide, in principle, information on the corresponding e-dependence for a rigorous relation to 
the MS factorisation scheme, a task which appears to be practically impossible. 

In order to circumvent this difficulty, we remark that R in Eq. (|19b|) is directly expressed 
as a function of the variable 7, and that the leading Pomeron singularity occurs because of 
the saddle point identification 7 = j L (a s (t) / u)) . It is then conceivable that such a singularity 
will be replaced by a much softer one if the effective anomalous dimension variable becomes 
7 — 1res(&s(t) /^)) at resummed level. A replacement similar to this one was used in anomalous 
dimension space in the study of the scheme-change of [10]. In our framework, we are able to 
ensure in general that 7 res is the relevant variable by assuming that the Qo MS normalisation 
change occurs in a fc-factorised form, i.e., by taking the "ansatz" 

,(MS)/^ _ [ ^7 -ft 1 f (Qo), 



2m 7_R(7,cj) 



2 The normalisation factor No takes NLa; corrections too, which however coincide [20] with those obtained 
by the known fluctuation expansion at e = 0. 
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where R is a properly chosen 7- and o;-dependent coefficient and /u> (7) denotes the uninte- 
grated gluon density in 7-space in the Qo-scheme. The latter is directly provided by the e = 
small-x BFKL equation, possibly of resummed (RGI) type. It is then clear that, at LLx level, 
R in Eq. takes the saddle point value R(yi,(a s (t)/u), 0), which therefore should coincide 
with the expression (|T9|) in order to reproduce Eq. $1}. Furthermore, the NLx expression |T4j) 
can be reproduced too, by a properly chosen O (a>) term in the expression of i?; and simi- 
larly for further subleading terms in the cu-expansion of R. In other words, Eq. (|2H|) can be 
made equivalent to Eq. (JUjl at any degree of accuracy in the logarithmic small- 2 hierarchy, but 
has the advantage that the effective anomalous dimension is order by order dictated by /i> , 
possibly in RGI resummed form, and is therefore much smoother than its LLx counterpart. 

By the argument above, we expect the u-expansion of the fc-factorized scheme-change (J23j) 
to be more convergent than (pPfj) in the small-x region. This encourages us to implement it at 
leading level (u — 0), in which we have 

/+00 
p(t-t')f^\t')dt' , (24) 
-00 

where 

^=L^^h' (25) 

is pictured in Fig. [T]a. For r > it is close to a 6-function and for negative r it oscillates 
with a damped amplitude for larger \t\ and with increasing frequency. The difference of g( MS ) 
and g(Q°) involves a weight function distributed around t' ~t which, convoluted with f^°\t'), 
includes automatically the RGI resummation effects of the Qo-scheme. A word of caution is 
however needed about the accuracy of (}2"4"|) in the finite-x region, where subleading terms in 
the w-expansion of R in Eq. are needed, and are left to future investigations. 
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Figure 1: (a) The function p(r) (solid) and its asymptotic estimates for r — > —00 (dashed). 
(b ) Sketch of the fastest convergence contour in the complex 7-plane of the integral in Eq. (j2*3j) 
for r <C — 1, including the cuts and two saddle points that provide the dominant contributions 
to p for very negative r. 



Even if A(r) = p(r) — 6(r) is in a sense a small quantity - - because the first two r- 
moments of A(r) must vanish — the numerical evaluation of (|24j) is delicate because of the 
large oscillations of p in the negative r region. For r — 1 the fastest convergence contour is 
shown in Fig. [TJb, where two saddle points 7 sp , 7* p of order % p (t) ~ 1 + iexp {\t\ + are 
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found. A saddle point evaluation 




<ro) + j; Q dr' 7sp (r') 
p(r) | sp = \l-lm{ 1 ° i > , (26) 

provides a good estimate for the contributions from the diagonal parts of the contour, while we 
integrate numerically the remaining part of the contour. The result (|2fi|) is To-independent, but 
the constant of integration c(tq) is determined numerically, e.g., c(— 3) ~ —1.18 + i 2.15. The 
function p(r) is also computed entirely numerically for r > —7. 

In order to illustrate the difference between small- x gluon distributions in the Qo and MS 
factorisation schemes, we consider a toy gluon density obtained by inserting a valence-like 
inhomogeneous term / in the RGI equation of [6], as follows 3 

f (x, t) = Ar°' 5 (l - x) 5 5{t - t ) (e'° = k 2 = 0.55 GeV 2 ) , (27) 

and solving the corresponding evolution for the unintegrated gluon density f(x, t) = G(Q 2 , k^; x) 
where t = log Q 2 /p 2 . The normalisation A is set so that the inhomogeneous term has a 
momentum sum- rule equal to 1/2. The solution of the RGI equation approximately maintains 
the sum-rule for the full resulting gluon, though not exactly because of some higher-twist 
violations. 

We then define the integrated densities 



xg (Qo \x, Q 2 ) = J dt' G(t - t') f(x, (28) 
xg^(x,Q 2 ) = J dt'p(t-t')f(x,t') , (29) 

which are shown in Fig. El in comparison to CTEQ (NLO) and MRST (NNLO) fits for the 
MS gluon at two scales. Note that we have not attempted to fine-tune the inhomogeneous 
gluon term to get good agreement at large x since in any case we neglect the quark part of the 
evolution which is likely to contribute non-negligibly there. 

One observes that the difference between the Qo and MS schemes is modest compared to 
that between the CTEQ and MRST fits and, in particular, there is no tendency for the MS 
density to go negative. This is despite the violently oscillatory nature of p, which might have 
been expected to lead to significant corrections. This can in part be understood from the 
approximate form of g( MS > 

xg^(x, Q 2 ) ee xg^(x, Q*)-^^ [xg™(x, Q 2 )] , (30) 

which is obtained by expanding the scheme-change in the anomalous dimension variable 7 ~ 
d/di to first non-trivial order. We can see that this approximation is pretty good in the small-x 
region for reasonable values of a s (Q 2 = 20 GeV 2 ), corresponding to an effective 7 ~ 0.25. 4 

We can also use our usual techniques [26] for extracting the splitting function itself in the 
MS-scheme. The results are shown in Fig. El where the MS curve is supposed to be reliable 



3 We adopt a variant of the NLLb resummation scheme introduced in Ref . [6] , where we perform the w-shift 
also on the higher-twist poles of the NLcc eigenvalue; we denote it NLLb' . The running coupling, as a function 
of the momentum transfer q 2 , is cutoff at q 2 — 1 GeV 2 . 

4 At the lower Q 2 value we do not show the MS3 approximation, because in general it breaks down for low 
Q 2 . 
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,0 



Figure 2: Toy gluon at scales Q 2 = 4 GeV 2 and 20 GeV 2 in the Qq and MS schemes, compared 
to MRST2004 (NNLO) [15] and CTEQ6M1 (NLO) [16]. At the higher Q 2 value we also show 
the MS3 approximation to the full MS evaluation. 




Figure 3: Resummed P gg splitting function in the Qq and MS schemes together with the 
uncertainty band for the Qq scheme that comes from varying the renormalisation scale for a s 
by a factor x^, in the range 1/2 < x M < 2; shown for two Q 2 values. 
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in the small- x region only (x < 10 -1 ), because at finite x the w-dependence of the scheme 
change will become important. It appears that the splitting function is more sensitive to the 
scheme-change than the density itself. At the lower Q 2 value the difference between the MS- 
scheme and the Qcr scneme (with NLLb' resummation) is nearly the same as the renormalisation 
scale uncertainty and so might not be considered a major effect. Recall however that the scale 
uncertainty is NNLx whereas the factorisation scheme change is a NLx effect, so that while the 
renormalisation scale uncertainty decreases quite rapidly as Q 2 is increased (right-hand plot), 
the effect of the factorisation scheme change remains non- negligible. 

At small x, most of the difference between the MS and Qo splitting functions seems to be due 
to the MS splitting function rising as a larger power of 1/x. One can check this interpretation 
by estimating the factorisation-scheme dependence of the asymptotic x~ Wc behaviour, using the 
MS3 approximation. We find 

Au c (t) = ^ C MS - ~ ^ ,y V ' y> ~ 0.029 , (31) 

dt 3 g'[t,u c [t)) 

where the numerical value is given for Q 2 = 20 GeV 2 , and is in rough agreement with what is 
seen in the left-hand plot of Fig. 01 



3 Approximate resummed splitting function for the MS 
quark 

Let us now discuss the inclusion of quarks in the resummed, small-x flavour singlet evolution. 
Resummation effects will be included via the unintegrated gluon density as discussed before, 
while the scheme-change to the MS-quark will be an (approximate) fc-factorised form of the 
one arising at first non-trivial LLx level. 

In order to better specify the scheme change, we first recall [18] that the quark density in a 
physical Qo-scheme corresponding to the measuring process p (e.g. F2 or Ft in DIS) — which 
we call p-scheme — is defined, at NLx level, by fc-factorisation of some g — > qq impact factor 
He P \ as follows: 

q®(t) = a»J d 2+2£ fc' ff E «(fc, k')T e {k') . (32) 
By then working out this convolution in terms of the eigenvalue function "K^ + e) 

(v) 

of H?', one directly obtaii 
dimension in the p-scheme: 



(p) 

of H?', one directly obtains a NLx resummation formula for the e-dependent qg anomalous 



[g^'Ht)]- 1 lgW(0 = 7 W(a s (i),o;; £ ) , (33) 

where, in the e = limit, 7^ = a B (t)JtQ (tl) has been obtained in closed form in various 
cases [18,20], including sometimes the e-dependence. 

The MS-scheme is then related to the p-scheme by the transformation 

qto = C$q m +C$g® iS > , (34) 

where we set Cqq = 1 and b = at NLx level, so that a s (t) = a A1 e e *. We thus obtain, by the 
definition (jHSJ) and by Eq. (JTJ), 

7 J)(a s ,,; £ )% )W ; £ ) = (^^[Cj)^,;^] +7 (f )(a B ,u;) 



lL (^)+eD] Cj 1 f(a s ^;e) + 1 ! 1 f s \a s ,u) , (35) 
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where £) = a s d/da s and 7<^ 1S ' ) is universal and e-independent, so that the process dependence 
of is carried by the perturbative, scheme-changing coefficient Cqg . 

The coefficient C$ in Eq. §^ can be formally eliminated by promoting e to be an operator 
in a s -space, as follows: 

e = - 1 ^)D-\ (36) 

v it) 

so that the square bracket in Eq. ()35|) in front of Cqg just vanishes (this procedure is rigorously 
justified in [20]). That implies that the l.h.s. of Eq. (JUS)) and ^F(t) are to be evaluated at 
values of e ~ 7l = C (a s /u) which modifies 7<^ 1S ' ) — 7^ at relative LLx level. Therefore, even 

if 7g!f S ^ is by itself a NLa; quantity, the subtraction of the coefficient part in Eq. (f33j) affects 
the scheme change at relative leading order, contrary to the gluon case of Eqs. (j!3|) and (|2*2*jl . 

The replacement (|3l)j) was used in [20] in order to get an "exact" expression of 7^ S ' ) at 
NLx level (that is, resumming the series of next-to-leading cu-singularities ~ a s (t)[a s (t)/u] n ). 
The main observation is that, by expanding 3-^(7) in the 7 variable around 7 = — e with 
coefficients Jin\e), and by converting Eq. (J32*j) to 7-space, we obtain 

±qjf\t) = a s (t) J d 7 6^(7)^(7,0;) 

f 00 , n \ _ 

^)+E^)^T fawfee^Mfru,)) , (37) 



n=l 



where ^(7,0;) is the Fourier transform of gi MS \t, 00) and "K{e) = 'Ni (— e) turns out to be the 
universal function 



T R 2 1 + £ 



_2tt 3 (l + 2e)(l + ge) 



e £ ^r 2 {l+e)T{l-e) 
r(l + 2e) 



J{ rat (e)J{ tran (£) , (38) 



which we factorise into parts with rational and transcendental coefficients. Note that the univer- 
sality of 9t(e), which is proportional to the residue of the characteristic function < K ( e\'y) / ^{'y+s) 
at the collinear pole 7 = — e, is due to the interesting fact that one can define [18] a univer- 
sal [20], off-shell g(k) — > q(l) splitting function in the collinear limit k 2 ~ I 2 <C Q 2 , for any 
ratio fc 2 /Z 2 of the corresponding virtualities. 

We then realise from Eq. (|37|) that the terms in the r.h.s. with at least one derivative d/dt 
are of coefficient type, and vanish by the replacement (|36|). The remaining one, proportional 
to 3i(e), is universal and, by (|3*6^1. yields the NLx result [20] 

7 « W( ),,)_, (i) ^(- 7L (^)^)g( 7L (^)^)^(^) , (39) 

where we have factorised a s (t) by the commutator [D, a s ] = a s and we have defined the quantity 

_ 00 _ 

n=0 

which has transcendental coefficients and differs from unity by terms of order (7l) 3 or higher, 
in the e ~ 7l region [18,20]. 

The complicated expression (f3T)|) has been evaluated iteratively [18], but not resummed 
in closed form. However, it can be drastically simplified in the fc-factorised framework by 
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neglecting the transcendental corrections O (7 3 ), on the ground that the resummed anomalous 
dimension is small. In fact, by setting R = J{ tran = 1 and 7l = a s /uj, Eq. reduces to the 
expression 

„, (i )r«f-^).l = ,(^f f^' ( (41) 



-(MS) 
<19 



u 1+D 



which is just the Borel transform of 



Tr 

47T 



1 1 1 

+ 



1 + 2e 3 1 + |e 



(42) 



Since IKJf = fa[2 n + 



1/2^ 
3^3' 



we obtain from (f4~T]) what we call the "rational" approximation 



NLx| rat (first derived in [18]) 



-(MS) 
>19 



OC B (t)T R ( 2 Q£ 1 2«S 



(43) 



This result can be further interpreted in fc-factorised form 



dt 



(MS) 
rat 



(44) 



by using the characteristic function 



a s T R 1 

47T 7 



e 2 ^ + -e§^ 



(45) 



which yields the result in Eq. at the LL saddle point. 

Finally, since the exponentials in (|45|) generate translations in t, our rough estimate of the 
resummed P qg is provided by the simple formula 



g(t + 2,x) + -g(t + -,x) 



(46a) 
(46b) 



By replacing in Eq. ()46a)l the resummed gluon density [6] and by performing the necessary 
deconvolution [26] we obtain the results in Fig. 131 compared to NLO and NLx [18] results. 

A proper comparison of the curves in Fig.|3]can be made only in the small- x region x < 10 _1 
because Eqs. fl3U|) and (J4T?)) do not include the finite- a; perturbative terms. We then notice that 

small-x resummation effects in q^ S ^ are sizeable even around x ~ 10 -3 and somewhat larger 
than the gluonic ones. They are anyway much smaller than the corresponding ones of the 
NLx result, thus showing that the resummed anomalous dimension variable is pretty small, as 
already noticed in the gluon case. We can also check how good the "rational" resummation is, 
by calculating the effect of the first O (7 s ) transcendental correction, which is also shown in 
Fig. It appears that the difference is indeed not large and anyway much smaller than the 
difference between the results of NLx| ra t in Eq. (}4*3j) and NLx [18], both shown in Fig. 0] 



To sum up, we have proposed here a fc-factorised form of the Qo — > MS scheme-change 
(Eqs. f!23l) and (|44jl ) which allows a convergent leading log hierarchy, because of the smoothness 
of the resummed anomalous dimension. Applying the leading scheme change to the gluon 
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Figure 4: The MS g — *> q splitting function for two values of a s and in various approximations: 
at two-loop [dashed], the RGI resummed rational one in Eq. ()46|) [solid] and with the addition 
of the first transcendental correction [dash-dotted], the rational NLx approximation [dotted] 
and the complete NLx one [dash-dot-dotted] 

density and — in an approximate way — to the quark density we have provided predictions 
for the gg and qg splitting functions in the MS-scheme and for the corresponding densities. 

We find that the gluon density itself is rather insensitive to the scheme change, while its 
splitting function is somewhat sensitive. Resummation effects for the MS quark are more 
important, but anyway much smaller than those at NLx level. We are thus confident that the 
scheme change can be calculated in a reliable way in a fully resummed approach as well. 
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